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!• IHOODUCTICN 



A. EACKGBCQKE 

As ct 1S74, the United States contained only 631 cf the 
world population but was using 33X of the energy consumed 
each year according to Ref. 1. Approximately one third of 
the oil used in the U.S. is imported. This heavy dependence 
cn foreign sources fcr energy coupled with the exhaustion of 
domestic fcssile fuels, particularly oil, in the foreseeable 
future is causing the D.S. military and civilian sectors of 
the economy tc search for energy sources that would be 
inexhaustible and independent of foreign control. 
Increasing fuel costs and environmental problems are other 
forces pushing the search for alternatives. Nuclear power, 
once considered the cure-all, is facing severe problems. 
Ccsts cf construction, operation, and fuel have risen 
dramatically. Questions concerning the safety of operation 
of the nuclear plants and of the storage cf nuclear wastes 
have become political issues that have slowed the 
exploitation cf nuclear power. The present nuclear reactors 
use a fuel that, like petroleum, has a limited availability 
and whose price is increasing rapidly. The breeder reactor, 
which would create more fuel than it burns, and the fusion 
reactor, which would use the hydrogen found in water, are 
still a Icng way from producing power for a world whose 
demands fcr energy increase each year. 

Hany people advocate the exploitation cf the sc-called 
"free" €nerc 7 sources. These sources are: (1) the kinetic 
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energy of moving fluids, such as winds and ocean currents; 
(2) the fctential energy of tides and rivers; (3) the heat 
generated inside the earth; (4) the direct conversion of 
solar energy into electricity and heat; and (5) the use of 
sun warmed water in conjunction with a source of odder 
water to provide the temperature differential to run a vapor 
power cycle. Sene cf these sources are truly ”free“. Shat 
is happening is the trading off of the transportation, 
processing, and environmental costs of a concentrated energy 
source for the capital, operating, and social costs of 
converting a very diffuse source into a mere concentrated 
ferm such as electricity. 

Ihe earth can he thought of as a giant heat engine, 
absorbing energy from the sun and reflecting energy as 
thermal radiation tack to space as pictured in Eef. 2. 
Since the pcles receive less energy than the equator, the 
atmosphere and ocean attempt to distribute the energy more 
evenly over the earth. The ait is heated hy absorption of 
sunlight and ty contact with the surface of the earth. The 
heated air rises thus setting up surface wind currents as 
ccclet air tries to replace the rising air. The shear 
stress created by the relative motion between the air and 
the surface water causes the water to move thus creating the 
surface water currents. Many currents in the oceans, 
regardless cf depth, originate with the wind shear stress at 
the air-surface interface. The circulation patterns cf the 
oceans beceme extremely complex due to the influence of the 
rotation of the earth and the shapes of the ccean basins. 

In general the surface waters near the equator are 
warmed hy the sun and flow towards the poles giving up 
energy by evaporation, and by convection tc the air and by 
radiation back to space. As these currents cool, seme of 
the water hecemes dense enough to sink and then flows back 
to the equator along the bottom. 
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The relaticnship between salinity and temperature and 
the density cf seawater is such as tc create a stable 
stratif icaticn of the ocean in many areas. A salinity 
decrease cr a temperature increase reduces the density. The 
total effect is to create a warm surface current overlaying 
the cold water returning from the poles. In the tropical 
waters the surface temperature may exceed 80°F f26.7oc) 
while 3,00C to 6,000 ft. (914 to 1828 m) dcwr the 
temperature nay be about 40®F -(4.400). Ibis temperature 
difference can be used to rur a man-made heat engine, like 
the Ocean Ihermal Energy Converion (OTEC) power plant 
ccEcept. 

I 

The proper choice of a site is critical to the cost of 
such a heat engine. Since the temperature difference is low 
compared to conventional thermal power generating methods 
which typically use temperature differentials exceeding 
SOOop (260OC) , the effect of a 1 ®f (0.56oc) loss in the AT 
is much more costly at the low thermal differentials, which 
are usually less than 50op (27.8oc) in the ccean. 
Therefore, the vertical temperature profile of the cceans 
will ccntrcl the location of an OTEC power plant. The 
surface temperatures vary with location and with time. At 
most places in the ocean there is a seasonal change of the 
surface water temperature of 5 to 8®F (2.8 to 4.4®C) except 
in the tropical areas that remain relatively constant 
according to Bef. 5. However, all regions of the oceans are 
affected ty wind-generated mixing of the surface water that 
changes the surface temperature by a degree or sc and a 
daily cycle cf a few tenths cf a degree. Some cf the 
possible 0.£. near-shore sites being considered are the 
waters around the Eawaian Islands and the Gulf Streai off 
the southeastern U.S. 

In the past, the engineering analysis and the cost 
estimation have teen conducted as if they were separate 
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fuECtions. In a system with a saall tempeiature 
differential, sqcb as this one, the engineering design and 
eccncmic analysis should be linked together as a single 
procedure because of the capital intensive nature cf the 
ocean thermal power plant. For example, the heat transfer 
egration (C = DA AT) is coupled to the ccst sguaticn for 
heat exchangers ($ = K a”) through the surface area of the 
tubes. Eecause of the small AT involved the area must be 
very large tc supply the same Q. The other components are 
all linked together by the engineering and ccst equations so 
that the problem cannot logically be broken into an 
equipment by equipment optimization. The present analysis 
considers a closed vapor power system operating cn the 
thermal potential available from the ocean. The economic 
optimization cf the system is carried out to show what 
information can be gotten, what conclusions can be reached 
from such analysis and where research effort should be 
expended tc improve the design of the system. This study 
does not attempt tc predict the cost of such a system or to 
make specific recommendations about the parameters. 

A schematic diagram of the system is shown in figure 1. 
The warm surface water, the Gulf Stream for example, is 
pumped thrcugh a heat exchanger. There some energy is 
transferred by boiling the working fluid such as ammonia. 
The working fluid is piped to a turbine where some of the 
energy is converted to mechanical energy then to electrical 
energy in a generator. The unavailable portion of the 
absorbed energy is rejected to the environment in the 
condenser and the working fluid is returned to the liquid 
state. The low temperature in the condenser is maintained by 
pumping cold ocean water up from the depths. Finally the 
working fluid is pumped back to the boiler. 

The Ccean Thermal Energy Conversion program has received 
much attention in the last several years from many 
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FIGURE 1. SYSTEM SCHEMATIC DIAGRAM 
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res€acch€rs funded at first primarily by the Haticnal 
Science Fcurdation (HSF) . Since the establishment cf the 
Energy Besearch and Development Administration (EBDA) , new 
research ccrtracts are funded mainly by this agency. The 
C.S. Navy is alsc funding some research in the OTEC area. 
Hef. 4, 5 and 6 give a good summary of the work that has 
been done. ihile much of the work has concentrated on the 
heat exchangers, the other components and the system as a 
whcle* are receiving attention. Some researchers are 
concentrating on the environmental and legal questions of 
such systems. Cther teams are considering hew to use the 
energy supplied from an OTEC power plant. 

In any design, trade-offs must be made. Capital cost 
can be traded for operational costs. Pumping power can be 
traded for heat exchanger size. Trade-offs can involve 
items that have a common measure such as dollars; seme other 
trade-offs are very hard to measure in common terms such as 
social and environmental costs. The trade-offs may involve 
crly one component such as a pump or a turbine, but most 
often the tiade-effs affect more than one component or 
possibly the entire system. 

The trade-eff process should continue until an optimum 
design is cbtained. from Ref. 7, "The aim cf optimization 
is the selection, cut of the multiplicity of potential 
sclutions, of that solution which is the best with respect 
to some well-defined criterion." In order for the plant 
design tc be optimum, the components must be optimized 
within the context of the entire system. Some cf the 
parameters fer a piece of equipment may have no direct 
effect cn the optimum design of another portion cf the 
system, tut seme parameters will have significant effects. 
For example, a heat exchanger is made up cf a large number 
of tubes cf a given diameter and length and is designed to 
have a certain heat transfer rate for a given flow rate of 
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water puaped through the tubes. The diameters, length and 
the number cf tubes may be optimized for lowest cost cf the 
heat exchanger. This approach leaves out the effect that 
the tube diameter and the tube length have cn the pressure 
drop across the heat exchanger. Minimizing the cost cf the 
heat exchanger involves minimizing the surface area. This 
may lead tc increasing the pressure drop across the 
exchanger. The higher the pressure drop the more powerful 
and the mere expensive the pump must be to maintain the flow 
velocity. Ibis is just one example cf the hundreds of 
trade-offs possible. 

Dnfortunately, what happens in many cases is that each 
group of experts working on a particular component optimize 
that item arcund some given conditions that may have been 
set by teams working cn some other section cf the sjstem. 
Little analysis is made of how the entire system responds to 
the coupling constraints. All that the pump manufacturer 
wants is that the buyer tell him what the flow rate, head 
and service conditions are and he will quote a price. 
Before the buyer can decide on the flow rate and head he 
should first know how the cost of the pump varies as flow 
rate and head vary so that those two parameters can be 
optimized in the context of the system. 

The pichlem then becomes one of acquiring sufficient 
infermatien cn hew costs vary as certain parameters change. 
Each manufacturer is able to make estimates cf costs fer his 
product tut he is reluctant to devulge his information to 
buyers cr competitors. This makes constructing cost curves 
a hazardous process at best, since the data must be gathered 
by other vetheds that may cause large unexplained variances. 

Sometimes the optimization analysis is made with 
assumptions that make the answers of questionable value. 
For example, sometimes the heat transfer coefficient, 0, is 
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assumed fixed in the equation g = OA AT. But, as shewn in 
Ref. 8 , for toiling the heat transfer coefficient is very 
sensitive to the temperature difference with the result that 
g = 0A'AI3.33, xhe accuracy of the engineering equations 
should he kept in mind. Most of the equations are 
correlations to fit experimantal data. The correlation was 
developed for some particular set of data and may disagree 
significantly when compared to data taken by other 
experimenters. The limitations of the various relations 
must be kept in mind when reviewing the results. For a 
system the size of the OTEC plant, a pilot plant should be 
constructed in order to prove out the answers given bj the 
calculations before expensive mistakes are made. 

As a design progresses, the analysis begins with a 
general picture containing many simplifying assumptions and 
proceeds foreward with more details added at each stage. 
This report is the second one written at the Naval 
Pcstraduate Bchccl on the subject of OTEC. Commander Furman 
Sheppard, DSN, in his analysis. Ref. 9, attempted to show 
what kinds of information could he gained from a combined 
thermal economic model. Be advocated the use of zone 
analysis, as developed by Hismer in Ref. 10, for large 
complex systems. In Hismer' s approach to zone analysis, the 
system is broken down into zones containing one or more 
cempenents. Ihe zones are connected by linking variables. 
Twe methods of optimization are possible. For the first 
method, in each zone, the zone parameters are optimized 
while the licking variables are held fixed. Then, the zone 
parameters are held fixed while the linking variables are 
optimized. Ihe process iterates between optimizing the zone 
parameters and the linking variables until convergence 
criteria are satisfied. In the second method the zones are 
cut off from each other and each linking variatle is 
separated into two variables, one on each side of the 
cutting place. The optimization proceeds until the linking 
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variables cn each side of the cuts are equal. 



Sheppard's system did not produce vapor in a boiler but, 
instead, maintained the working fluid in the liquid state in 
the warm seawater heat exchanger and vaporized the working 
fluid, attJicnia in his model, by expanding it through a 
“black hex" called a vaporizer. His optimization analysis 
was confined to the zone containing the warm seawater 
exchanger, the feed pump, and the warm water circulating 
pump. 



E. CBJECIIVIS 



The present analysis improves the previous model 
analysis hy replacing the liquid-to-ligui d heat exchanger 
and the vaporizer with a boiler. The boiler model is to be 
of sufficiert detail so that a realistic design is 
represented. The ether major components of the system, such 
as the turbine, the cendenser, and the condenser circulating 
pumps, are included, k nonlinear programming technique is 
used to perform the thermal-economic optimization of the 
completed system. The research reported here is intended to 
show the usefulness of optimization analysis in the design 
of a low thermal difference system and to demonstrate some 
of the problems that must be considered in this type of 
analysis . 
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II. THE MODE L 



A. 6ENEB1L EESCBIETION 



Vapcr pc«€x generating cycles usually consist of four 
processes. Energy is absorbed from a high temperature 

I 

source and is used to vaporize a working fluid in a boiler. 
Next, work is extracted by expanding the vapcr in an engine, 
exhausting tie vapcr at a lower temperature and pressure. 
Th-e working fluid is then condensed back to the liguid state 
by rejecting energy to the Icwer temperature energy sick in 
the third process. Finally, the liquid is returned to the 
higher pressure of the boiler by a pump. In the following 
analysis, the energy source is a current of ocean water 
heated by the sun and the energy sink is assumed, to be 
colder ocean water from deep currents. The working fluid is 
assumed to be propane although many other fluids are 
possible, such as ammonia or the Freons. 

Six ccmpcnents are modeled in the follcwing analysis. 
The hot side heat exchangers, the boilers, vaporize the 
working fluid which is expanded through a turbine exhausting 
to a condenser that rejects the unused heat to the cold 
sink. Three pumps are included, a feed punp returns the 
working fluid to the boiler and two pumps circulate the warm 
and cold waters through the boiler and the condenser, 
respectively. The arrangement of the system is shown in Fig. 
2. The nodes are numbered to correspond with the fluid 
state points shown later, in Fig. 3. 



21 



CJ 

H- 1 

H 

< 




22 



FIGURE 2. SYSTEM ARRANGEMENT 



of the 



The analjtical model consists of the equations 
describing the system and a set cf parameters to be 
optimized sc as tc minimize the capital cost of the plant 
subject tc certain constraints. The objective function, the 
function tc be minimized, is the sum of the capital costs of 
the pumps, the boiler and the condenser. Although the 
turbine is a significant capital cost and mould affect the 
design cf the system, it is not included in this model. 
There are two classes of constraints. The first class is the 
explicit constraints that specify the upper and lower bounds 
on the parameters. The other class of constraints is the 
implicit constraints. These are the engineering equations 
that describe the operation of the power plant. For 
example, a ccnstraint might say that the temperature inside 
the boiler must be less than the temperature of the warm 
seawater because the equations fail if the requirement is 
not satisfied. Another example would be a power balance 
that requires that the sum of the energy flows into and out 
of the system be zero. Additionally, they proscribe the 
region that the optimization routine may search within for 
the minimum cf the problem. The implicit constraints are of 
two types, equality and inequality constraints. Inequality 
constraints set one-sided bounds on the feasible region. 
Some of the inequality constraints are restrictions placed 
on the maximum or minimum values of the parameters due to 
judgements made cn the basis of practical engineering 
considerations. The equality constraints are relationships 
that must be strictly satisfied at the optimimum. The system 
of equations describing the model fit into the following 
form: 
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f (X) 



subject to: g (X) > 0, i=l,2,....,m 
i 

h (X) = 0 , i=ni+ 1 , . . . . ,m+n 
i 

wheie X is tfce vector of the parameters tc be optimized, 
f(X) is the capital cost function, g^ (X) are the inequality 
ccnstraints and hj (X) are the equality constraints. The 
parameters are the dimensions of the toilers and the 
condensers, the boiling and condensing pressures, the sea 
water velocities, and the mass flew rate of the working 
fluid through the system. 

The tasis for the economic framework must be set cut at 
the start. ill measurable costs ace relevant to the final 
design of the power plant, but, for some types of analysis 
and design, many costs can be ignored. The life cycle cost 
of the CTEC must be considered before the decision is made 
to commit the O.S. and the rest of the world tc the 
consequences cf chcosing this energy generation scheme as 
even a partial solution to the energy problem. The life 
cycle cost consists of the capital costs, the financing 
costs, and the operating costs. The trade-eff of operating 
cost and capital cost is not considered in this research, 
because the operating cost is, for the most part, 
independent cf the overalL dimensional characteristics of 
the components. Operating costs are deminated by the 
structure, site location, and maintenance requirements. 
Such things as the kind of instrumentation, controls, 
bearings, and structural materials have mere effect cn the 
operational ccsts than the length of the boiler or the 
diameter cf the toiler shell. Therefore, cnly the capital 
ccsts of the major components, except the turbine, are 
considered. 
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To b€ cseful for model analysis cost data lust be 
transforiBed into a cost estimating relationship (CEB) that 
shews hew cost varies with some variable or some combination 
of variables. The costs must be demonstrated to be highly 
correlated with the variable chosen. Since the cost 
relation is cnly a correlation, the bounds on the accuracy 
of the eguaticL should be specified by those who develop it. 
Since the relation is arrived at from a finite body of data, 
the CEB is directly applicable only in the range of the data 
collected. Eitrapolation of the CEB outside the range is 
dangerous because there is no data in that region to support 
the assumed curve. However, the value of CEB's is the use 
of them tc enable the estimator to make a better analysis of 
a new design than he could if the information is not 
available. The relevant range of the CEB must be kept in 
mind wher analyzing the solutions from the standpoint of the 
reliability of the answers. Another factor to consider when 
using cost estimating relationships is whether the equipment 
the data was taken from is comparable to the system the 
analyst is considering. If the design, the technology, or 
the application of the proposed equipment is significantly 
different from that from which- the data is taken, the CEB 
may give unreliable answers. 

This analysis begins with the development of the 
engineering equations that constrain the system. This is 
followed by the development of a cost framework based on 
capital costs as the function to be minimized. The 
engineering and economic relations are linked together in 
the objective function and the constraints. A description 
of the algorithm that performs the optimization is included. 



B. CYCLE AHAIYSIS 



The tuilcing of the model begins with the analysis of 
the ther ncdyna,mic cycle. Figure 3 is a t emf erature-ent ropy 
diagram cf the working fluid on which the cycle is shown. 
These nuitbered nodes represent the thermodynamic state 
points at that location in the system and correspond tc the 
numbers ic Fig. 2. The boiling and condensing processes are 
assumed tc be reversible processes. In the pump anc the 
turbine the reversible process is indicated by the subscript 
"s" on the state pcint number. In the actual system there 
are irreversible losses inside each component that are left 
out of the analytical model or lumped in with ether losses. 
The pressure drops due to frictional losses in the toiler 
and the cendenser are neglected. 

State pcint 1 assumes that the liquid leaving the toiler 
is saturated. State points 2s and 2 represent the 
conditions leaving the feed pump. If there were no losses in 
the pump, the fluid would, be at 2s. ‘ The fluid enters the 
bciler tc te heated, first, to saturation conditions at 
state point 3 then boiled at essentially constant pressure. 
State pcint 4 shows the vapor leaving the bciler as a 
saturated vapor. The lines from point 4 to 5s and 5 
represent the expansion through the turbine. The path tc 
state pcint 5 accounts for the losses due tc friction and 
other irreversible processes in the turbine. The remaining 
heat is then rejected from pcint 5 to point 1, again assumed 
to te a constant pressure process. 

The theoretical heat transfer rate in the boiler and the 
cendenser is found by subtracting the enthalpy at state 
points 2 and 1 from the enthalpy of state points 4 and 5, 
respectively, and multiplying by the total mass flow rate of 
the working fluid. 



Qab - - Ha) 
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FIGURE 3. THERMODYNAMIC CYCLE 
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Qrcj = Mw^(H5 - «1) 



( 2 ) 



Ih€ enthalpy at state point 5 is determined by first finding 
the moisture content of the fluid exiting the turbine. 



The enthalpy at pcint 5s is determined from 



<3) 



fls, = + i 



This allc«s fig to be found from 






- H5 



(h) 



t5) 



where the efficiency of the turbine is assumed kncwc. The 
power out of the turbine is now known. 



Ptur* = - Hj) (6) 

For the feed pump the change in enthalpy is accounted for by 
the change ir pressure assuming negligible rise in the 
temperature. The enthalpy of point 2 is determined from 



. - Pi I 

- 0, 

where the efficiency is assumed known. The power required to 
drive the puop is determined from 
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- H^) 



( 8 ) 



Tc gain a tetter feeling fcr the low tbermal potential 
systeiES and tc icake comparisons with other systems several 
efficiencies are worthwhile calculating. Ihe first is the 
Carnot efficiency which defines the maximum efficiency 
pcssihle for a system operating between two temperature 
liiits in degrees absolute, (®K). 



fear 



= 1 - 



"^cold 

^hot 



(9) 



The next efficiency calculation takes into account the 
theoretical thermodynamic cycle and is the ratio- of the net 
poser out of the cycle to the beat taken up by the cycle 



~ , 10 ) 
QqIs 

Ibis acccunts for the thermodynamic inefficiencies inherent 
in the cycle. Finally, the system efficiency takes into 
account the other internal consumers of power: the 
circulation pumps, the pipe losses, and any other losses, 
such as the antifouling equipment. This calculation 
acccunts fcr all of the mechanical and electrical 
ineff iciercies of the system. 




NET POMES AT THE EUS EAE 
Q ob 



( 11 ) 



For example, if the hot water is at 75°F (23.9°C) and 
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the cold satec is assumed at 40®P (4.4°C) the Carnot 
efficiencj is 6.53. For one typical set of assumptions, as 
shewn in Tafcle 5, where the turbine and the pumps are 
assumed to he 853 efficient, the cycle efficiency is 2.63 
and the system efficiency is 2.43 for one feasible design. 



C. ENGIMZiBIHG FRAMEWORK 



This section describes the design of the components that 
make up the system. The level of detail of the analysis 
provides a physically realistic model and closes the loop. 
Ite pipe friction losses on the cold and hot water inlets as 
well as the pipe losses in the working fluid side are 
neglected. The turbine is included only as an energy 
transducer and no attempt is made to obtain costs, physical 
dimensions, cr fluid flow characteristics. The efficiencies 
of the turbine and the pumps are assumed as known constants. 
The heat absorbed in heating the working fluid to the 
saturated state is a small portion of the total energy taken 
up in the boiler; therefore, the sensible heating cf the 
fluid is neglected. 



^ Ec iler 



The bciler consists cf a number cf tube-in-shell 
beat exchangers connected in parallel as shown in Fic. 2. 
Figure 4 shows the internal arrangement of one shell. The 
tube bundle is horizontal and boiling occurs on the outside 
surface cf the tubes. The tube bundle is submerged in the 
working fluid. The working fluid is maintained at a level to 
just cover the top row of tubes. The warm seawater is pumped 
through the inside of the tubes. 
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FIGURE 4. BOILER CROSS SECTION 
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Ih€ clearance between the tubes is fixed at 0.5 
inches (1.27 cm) to allow for flaring of the tube ends or 
ether space and strength reguirements. Similarly, the tube 
thickness is fixed at 0.028 inches (0.071 cm), since no 
strength reguirements due to differences between internal 
and external pressures on the tubes or material loss due to 
corrosion of the tube surfaces are accounted for. 

At the tep of the tube bank the leaving vapor 
contains a large amount of liguid droplets which must be 
eliminated either in the boiler or by an external separator. 
In the mcdel the separation takes place in the boiler by 
allowing sufficient space above the tube bank. Falen and 
Small, Bef. 11, presented a method for finding the necessary 
separation volume that considers the effects of surface 
tension, density of the vapor and the liguid, and the mass 
flew of the working fluid. First the allowable vapor lead is 
found from 



VL = 2290A f 1°*^ (12) 

16.86x10-5 )J 

Taper load is an expression for the mass flow rate of vapor 
through a unit volume. The expression for the maximum vapor 
lead limits the flew rate so that the moisture droplets are 
able to settle out of the vapor before it leaves the bciler. 
Then the reguired segmental area is found from 



SA 



Mu,f 

VL L S 



(13) 



The segmental area is the cross sectional area from the top 
of the tube bank to the top of the shell. Knowing the 
segmental area the height of the tube bank above the main 
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diameter is found from 



2 

SA = C,25 Cs ABCOS 




2 2 0. 5 

hs(0.25 Es - bs ) (14) 



Palen and Small also present a much simpler method of 
determining the total height of the tube bank. The height 
of the separation volume is set at 40X of shell diameter 
whioh results in "hs” being equal to 0,1 Ds, Either method 
can be used alone or the first method can be used as a 
constraint on the second. For this problem the segmental 
area is found by both methods uith the first used as a check 
on the values given by the second. 



tube, the number of tubes per shell must be determined in 
order to find the total energy absorbed by the heat 
exchanger. The tubes are assumed to be laid out in a 60 
degree triangular pattern for maximum compactness. The 
tubes are located at the vertices of an equilateral 
triangle. Two cross sectional areas are calculated. The 

area of a hexagon whose minor radius is equal to one half 
tbe pitch is determined from 



Then the effective cross sectional area of the tube bank is 
found by subtracting the clearance area, SA, from the cross 
sectional area of the shell after correcting the shell 
diameter for a clearance space to allow for free flow 
channels both within the tube bundle and around the outside 
of the tuhe bundle. The effective cross sectional area of 
the tube bank is divided by the area of the hexagon, found 



Since the beat flux calculation is based on a single 



A 




(15) 
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afccve, tc determine the number of tubes in a shell. 

Ite flew of the liquid and vapor through the tube 
back is assumed to be turbulent. The incoming subcccled 
working fluid from the feed pump mixes with the 
recirculating fluid returning frem higher up in the tube 
bank. The circulation of the fluid is induced by the 
boiling process. This mixture enters the tube bank and 
begins the boiling process. Near the bettom of the tube 
back the teat transfer goes mostly to bring the fluid up to 
the boiling temperature. After the fluid reaches saturated 
liquid cenditiens, stable boiling begins. The vapor bubbles 
rising threugh the surrounding liquid increase the 
turbulence around the tubes. As the liquid continues up 
threugh the tube bank, the local boiling temperature is 
decreasing due to the decreasing hydrostatic head. This 
decrease in the boiling temperature increases the 
temperature difference and increases the heat flew and 
increases the bubble generation rate. The vapor generation 
may become sc rapid that a condition called vapor blanketing 
may ccccr in the interior of the tube bank where the liquid 
cannot flow inward fast enough to displace the vapor flowing 
through the region. This is similar to burnout in 
ccnvective bciling in that the tubes dry out and the heat 
transfer rate decreases. This effect becomes pronounced 
cnly near the critical heat flax according to Starezewski, 
in Bef. £. 

fill cf the heat and mass flow variations in the tube 
bundle are averaged by performing the calculations for a 
single tube located at the mid height of the tube bank. 
This tube is assumed to represent the average of all the 
tubes in the tube tank. The boiling situation is assumed to 
be saturated pool boiling. Quoting frem Collier, Ref 12, 
"pool boiling is defined as boiling from a heated surface 
submerged in a large volume of stagnant liquid". Since the 



34 



fluid vclume seen by the tube is neither large nor stagnant, 
a ccrrecticn is needed for the effect cf being in a tube 
bundle. Saturated boiling is defined as boiling at a 
constant temperature and pressure as the liquid changes to a 
vapor. Ihis is probably not strictly true because the 
boiling temperature is controlled by the local pressure at 
each tube, and the fluid is not reaching a stable 
temperature as it rises in the bundle. The situation is 
prcbably closer to convective boiling since the fluid is 
boiling in a confined channel between the tubes. An 
additiocal effect in the boiler is the fact that the 
hydrostatic bead at the bottom of the tube bundle causes the 
boiling tc taXe place at a higher pressure at the bottom 
than the tciling at the top of the bundle where the pressure 
is equal to the boiler exit pressure. Ihis effect is 
accounted for by using the average head of the bundle in the 
heat transfer equations. 

Ihe heat transfer analysis of the boiler begins with 
the basic equation 



Q = OA AT^ (16) 

where U is the overall heat transfer coefficient, A is the 
area through which the heat is transferred and is a 
properly defined mean temperature difference. In this 
paper, the product (JA is treated as a single quantity called 
the thermal conductance. One may look at the equation as an 
analogue cf the basis electrical equation I = E/R where I is 
equivalent tc Q, E is the same as AT^ and UA is 1/S. Ihis 
analogy is used when developing the expression for UA. 

In order to account for the effects of the tube 
seeing many nearby neighbors, the equation is modified by 
the inclusion of a factor, F, to account for the effects of 



35 



tb€ surrcundicg tubes. Equation 16 is changed to 



Q = OA F AT^ (17) 

where F is defined as the ratio of the heat flux of the tube 
hurdle tc the heat flux for a single tube in a stagnant pool 
boiling situation. According to Sef. 13, a report published 
ty workers at Heat Transfer Research, Inc. (HTEI) of 
Alhambra, California, F is greater than cne. In Figure 5, 
which is sketched from a figure in their report, F ranges 
frcB 4.73 at AT„, = 8°F to 2,63 at AT^ = 40®F. Since nc data 
is given cn the conditions of the test, F is assumed to be 
fixed at three. 

lE order to find AT»»„ the heat exchanger 
effectiveness, 6 , is determined. The effectiveness is 
defined as the ratio of the actual heat transfer rate to the 
maximum possible heat transfer rate as shown by 



- ~ '^h,out ^ 

•“min “ ^c,in ^ 



where 



= minimum of or C^. (20) 

In the boiler, T^ is the tempera-cure of the entering 
seawater, acd T^/m is the saturation temperature of the 

working fluid at the pressure at the mid height of the 
boiler. Reference 14, by Kays and London, contains 

equations for finding the effectiveness fcr various flow 
configurations. In addition to finding must be 
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LOG q (BTU/hr-£t 




VERSUS A TUBE IN A BUNDLE 
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For the 



kncwn in crder to decide which equation applies, 
case of hciling, is infinite; therefore, regardless of 

the flew pattern 



£ = 1 - EXP(-NIO) (21) 

where NTD, the number of heat transfer units cf a heat 
exchanger, is defined as 

NTO = Ua/C,„in (22) 



and 



Cmin ~ 



(23) 



AT^is found frem 



OA AT,y. 



(24) 



and the temperature of the sea water at the outlet cf the 
tube is determined frem 

AT^ 



where 



AT« = T 



h,in 



- I 



Cjin 



(26) 



Ihe quantity UA is determined next. Figure 6 
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FIGURE 6. THERMAL NETWORK FOR A 
SINGLE TUBE 
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illustrates the uethod used. If OA is pictured as the 
reciprocal ci resistance then the flow cf energy froa one 
fluid tc the ether can be pictured as a series electrical 
circuit and OA is determined from 

OA = V(R1 + E2 + S3 + E4) (27) 

Beference 15, by Holman, shows how each resistance can be 
feund. HI is the thermal resistance due tc convecticn on 
the inside cf the tubes, the sea water side. 



B1 



1 

ii 1 



B2 is the thermal resistance of the tube wall. 



(28) 



In (do/dj ) 
2 7T L 



(29) 



B3 is the thermal resistance attributed tc boiling cn the 
outside cf the tubes. 



B3 



TTh^^ d, L 



(30) 



B4 is the theraal resistance due to any other resistances, 
such as feuling and corrosion deposits. 

Ihe heat transfer coefficient on the seawater side, 
h^^ , is determined from the dimensionless group called the 
Nusselt number defined as follows for flow inside of tubes. 
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Nu 



(31) 




In order tc lind the Nusselt numter, three distinct regions 
of flow are allowed for in the pzcfclem; laninar, 
transiticE, and turbulent flow. The division points are in 
terms of values of the Reynolds number. Re = ud/J/yu , 

laminar: 0 < Re < 2,0C0 

transition: 2,000 < Re < 10,000 

turbulent: 10,000 < Re < infinity 

In the laminar region the Sieder-Tate correlation is used. 



1/3 

Nu^ = 1.86 (Re Pr dj /L ) (32) 

In the turbulent region the Cittus-Boelter correlation is 
used. 



0.8 0.3 

Ku^y = 0.023 Re Pr (33) 

In the transition region the flow is unstable, but from data 
in Ref. 14 it appears that a fairly smooth curve drawn from 
the lamirar curve to the turbulent curve adaguately 
represents the value of the Nusselt number in the transition 
region. This is done by assuming the form 

= y Nu^ + (1 - y) Nuy^ (34) 



where 



3 2 

y=aRe+bRe+cRe+d 



( 35 ) 
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Th€ coeffici€Ets in the polynomial are found by requiring 

Nuj^ = Nu^_ at Re = 2,000 
Nu-p^ = at Be = 10,000 

and that the first derivatives with respect to the Reynolds 
number match at Re equal to 2,000 and 10,000. 

Ihe I. I. Mostinski correlation, from Ref. 16, for 
single hcrizcntal tube pool boiling is used to find the heat 
transfer coefficient on the outside of the tube. BTEI’s 
extensive experimentation led Palen, larden, and Tabcrek, in 
Ref. 13, tc conclude that Mostinski’s equation produced the 
•'acst consistent results”. His correlation is 



C.69 0.7 0.17 1.2 1C 

h^^ = 0.CC658 p^ q (1.8 r + 4 r + 10 r ) (36) 

where 

r = F/Pc 

= critical pressure of the working fluid (psia) 
g = heat flux (BT0/hr-ft2) 



A final heat flux relation needed is the relation fcr the 
maximum beat flux allowed. Hostinski,s correlation fcr a 
single tube is used. 



0.35 0.9 

= 003 p^ r (1 ” r) (37) 

where the units are the same as in equation (36) . From 
HT5I*s data, under unknown conditions, q„,<,x bundle 

appears tc be cne third of the ^ single tube. 
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Ite power required to pump the salt water through 
the tubes depends on the pressure drop across the tube bank. 
There are three causes of the pressure drop. The flow 
experiences a loss as the water enters into and exits from 
the tubes in addition to the frictional loss along the tube. 
From Bef . 14, the pressure drop fcr a liquid is determined 
from 



2 4 f L 

Ap = 0.5 Uj^ ( + Kc - Ke) (38) 

In order to determine f, Ke and Kc, the flow regime must be 
specified. The regimes are a little different than those for 
the Nusselt rumber. 

laminar: 0 < Re < 2,000 

transition: 2,000 < Re < 5,000 

turbulent: 5,000 < Re < infinity 

The Fanning friction factor, f, in the laminar regime is 

f = 16/Re (39) 

In the turbulent regime, the Fanning friction factor is 
found from the Blasius equation in the range 5,000 < Be < 
100,000 



0.25 

f = 0.079/Re (40) 

which is a close approximation to the Karman - Nikcradse 
equation , 
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-0.5 0.5 

(4 f) = -0.8 + 0.87 ln(R€(4 f) ) 



< 41 ) 



which is valid thrcughout the turbulent regicn (Re > 5,000), 
assuming smccth walled tubes. Either equation (40) or (41) 
may be used with similar results but the Elasius is tetter 
suited to this analysis because the Karman-Nikuradse 
equation mcst be solved by iteration. Equation (41) 
consumes more computer time to solve and can cause accuracy 
problems if not enough significant digits are evaluated. 
The fricticr factor in the transition region is determined 
by fitting a least squares curve through a hand drawn curve 
frcm the end of the laminar curve to the beginning of the 
turbulent curve. 

The entrance effect, Kc, and the exit effect, Ke, 
depend on the Reynolds number, the free flow to frontal area 
ratio, 4> , and, in the laminar region, on the tube length. 
For the boiler, is determined by 






2 

d,' H 



2 

Ds - 4 SA 



(42) 



Figure S, page 7-13, in Ref. 17 is used to find Ke and Kc. 
The effect of Reynolds number and L/d ratio was eliminated 
frcm the relationships in the laminar and the turbulent 
regimes of the flow, tut the Reynolds number was included in 
the transition region because of the large variation of Kc 
with Reynolds number. The following equations are used to 
find the entrance effect in the various flow regimes. 

laminar: (using the curve for 4 (L/d) Re = 0.1) 



Kc,_ = 1. - 0.4 



(43) 
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m 



lurbulent: (using the curve for Re = 5,00C) 



= 0.52 - 0.4 (j) (44) 

Iransiticr : 

= A KCl + (1 - A) Kc-p^ (45) 

where 

A = e Re + f (46) 

The coefficients e and f are determined tj reguirirg the 
transition Kc to egual the laminar Kc at Re egual to 2,000 
and the turiulent Kc at fie egual to 5,000. Ihe exit effect 
is the saae for all flow regimes, and is found frcn the 
following eguation that approximates the 4 (L/d) Be =0.05 
curve in Bef. 17. 

Ke = 1.C093 - 2.5178(^) + 1.1613 <}> - 0. 17677 <j) (47) 

2* Ihe Conde n se r 



The design of the condenser is almost the sane as 
the design of the boiler. It is a fixed tube sheet, 
tube-in-siell heat exchanger. The condenser does net need a 
separation volume so the entire shell is filled with tubes. 
In order to allow the entering vapor to reaoh all of the 
tubes, ohannels through the tube hank are needed. The 
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allowance for this is accomplished by decreasing the shell 
dianeter by a fixed factor when calculating the numter of 
tubes. The calculation of the working fluid side heat 
transfer coefficient, h , is different. The equation 
accounts for the condensate from one tube falling on those 
below it atd for the subcooling of the condensate. a.M. 
Chen's equation, as developed in Bef. 18 and shown in a more 
siiple form in Bef. 17, is used. 

0.25 

(48) 



= C.726 



/ J L n d, AT J 



where 



h = BlO/ftz-hr-op 



~ *^Sat ~ 





. 1 CpAI 



n = number of tubes in a column 



(4S) 

(50) 



H* is the enthalpy corrected for subccoling of the 

fg 

ccrdensate. Equation (48) replaces equation (36) it the 
analysis of the condenser. In the circular shaped tube 
bark, n is fcurd by determining the number of tubes in a 
column of average height. 



Ds 

n = 9.425 — - 1 (51) 

pi 



The definition of AT^ becomes 
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* 







(52) 



i h 

in contrast tc equation (26) above. The seawater cutlet 

temperature is found from 

which replaces equation (25) . 

In crder tc find the heat transfer coefficient, the 
temperature cf the outside tube surface is needed in 
equation (49) ; therefore, the thermal network is solved for 
the outside wall temperature and for the heat flow rate. 

a = Di> (T„ - ) (511) 



where 



OA* = 1/(R1 ♦ R2 + R4) (55) 

^bulk “ averaqe of the inlet and 
the outlet seawater temp. 

There is nc factor, F, in the heat transfer equation for 
condensation since the effect of neighboring tubes is 
accounted fcr in equation (47) . 



3 • Tte_gumpg 



There are three pumps in the model. Both 
circulaticn pumps for the boiler and the condenser are 
assumed to he propeller pumps (axial flow pumps) . The feed 
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puiBfs are assuned to be centrifugal pumps. Ihe differences 
in the characteristics of the two types of pump show up only 
in the calculation of capital costs of the pumps ir the 
model. 



Ihe choice of the type of pump to apply is made on 
the basis cf the flow rate, GFH, and on the head the pump is 
working against. Cualitatively, the propeller pump is used 
where the flew rates are high and the head is low. Ihe 
centrifugal pumps apply where the flow rate is comparatively 
lew and the head pressure is high. The gualitative picture 
is placed in tetter perspective through the calculation of 
the specific speed from 



Nsg 



0.5 

BPH G£H 
0.75 

Head 



(56) 



The centrifugal pump is used when 500 < Nsg < 7,500 and the 
propeller puop applies when 7,500 < Nsg < 15,000, according 
to Bef. IS. The choice of the type of pump is checked after 
the optimi 2 aticn process has been completed by solving 
equation (55) . If the value of the specific speed is 
outside cf the range, 500 < Nsg < 15,000, the maximum 
efficiency cf the pump drops off rapidly; therefore, few 
pumps are designed to operate outside cf this range. 

Since the pumps are "black boxes", only the power 
consumed by the pumps is calculated. The power required to 
drive the puop is found from 



P 



M Ap 
*?pum p 



(57) 
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c. 



CCSi: lEAHEROEK 



Is discussed at the start of this report, the capital 
costs of tie major components are the only costs to be 
considered ir this mcdel. The boiler and the condenser are 
both tu te-in-shell heat exchangers assumed to be cf the 
fixed tube sheet type. The only difference being that the 
shell of the toiler is partly filled with tabes. Both 
seawater circulation pumps are assumed to be propeller type 
, puaps and iiade cf the same material. The size and number 
per heat exchanger shell may be different fcr the bciler and 
the condenser. The number of working fluid feed pumps, 
which are the centrifugal type pumps, are assumed to equal 
the number cf bciler shells. 

The cost estimating relationships fcr the heat 
exchangers and the centrifugal pumps are taken from Eef. 20, 
by K.M. Guthrie. The eguation fcr the cost of the heat 
exchangers is 



Jhe = 11S-3 (F,|* Fp) F„ |58) 

where the ccrstants are calculated from a leg-log plot of 
cost versus total tube surface area and 

=0.6 for a fixed tube sheet 
Fp = O.C for pressure correction (p < 150 psi) 

F^ = material factor from Table 1 
Ihe ■ cest index 
= do L N 

fcr a relevant range cf 100 < A < 10,000 ft^. 

The cest eguaticn for the centrifugal pumps and the 



^9 



asscciated meters is 



0.721 

fm ?o I+i, Ch (59) 

where the constants are found from a straight line 

approximaticD cf the log-log plot cf cost versus C for the 

section cf the curve from C = 30,000 to C = 300,000 for an 
electric motor driven pump and 

Fq = 1.0 for suction temperatures < 250®F (121°C) 

= cost index 

F,y, = material factor from Table 2 
= prcdcct of the flow rate in GPM 

and the pressure differential in psi 

for a relevant range of 30,000 < < 300,000. The cost 

data in Guthrie's article is given fer a time base of 

isid-1968 . 

Reference 21, also by Guthrie, gives cost data for 
propeller pumps and centrifugal pumps. The eguation for the 
propeller pumps with the motors is 



0.783 

= 5.767 F^ Ff, Flow <60) 

where the cerstants are found from the log-leg plot cf cost 
versus flew rate for an electric motor driven pump and 

Fp = 1.0 fer a suction pressure < 150 psi 

F^ = material factor from Table 3 

Icp “ cost index 

Flow = flew rate in GPM 

for a relevant range of 1,000 < Flow < 100,000 GPM. The 
time base fer the propeller pump data is fer the erd of 
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Table 1 

HAIIEIAI FACTORS FOB THE HEAT EXCHANGERS 



M ateria l 



Shell 


Tube 




Fm 


Carbcn Steel 


Carbcn Steel 


1.00 


Carhcn Steel 


Brass 




1.52 


Carbcn Steel 


Stainless 


Steal 


3. 52 


Stainless steel 


Stainless 


Steel 


4.50 


Carbcn Steel 


Monel 




3. 75 


Mcnel 


Monel 




4.95 


Carbcn Steel 


Titanium 




11.10 


litatiuffl 


Titanium 




16.60 



Table 2 

BATEFIAL FACTORS FOR THE CENTRIFUGAL FEED PUKE 



Material 


Fm 


Cast Iron 


1.00 


Bronze 


1.28 


Cast steel 


1.32 


Stainless Steel 


1.93 


Monel 


3.23 


Hastelloy C 


2.89 


Titanium 


8.98 
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Table 3 

HATEBIAL FACTORS FOR PROPELLEE PUMPS 

pumps 



Material 


Fm 


Cast Iron 


1 .00 


Cast Steel 


1,28 


Stainless Steel 


1 .64 



1S70. Ih€ costs arrived at from equations (58), (59), and 
(6C) ace in dcllars per component, for example, dollars per 
heat exchanqer shell. The purpose of the cost index is 
allow all the equipment to be costed on a common basis. The 
time difference is assumed to be small enough so that the 
cost index is set to 1.00 for all the components. 

E. ASSEKELl Cf THE MATHEMATICAL MODEL 



The enqineering and the cost relationships are linked 
together in the model. The cost equation is the function to 
be mininixed and the implicit and explicit constraints are 
the engineering equations that restrict the feasible region 
of the prctlem. The explicit constraints restrict the 
variables, the vector X, to positive values. The rest of 
the model is formulated in the following form 

minimize f (X) 

subject to: g (X) > 0, i=1,2,....,m 

i 

h (X) = 0, i=m+1 , . . . . ,m+n 
i 

where the definitions of the elements of X are shown in 
Table 4. Several of the variables must have only integer 
values if the system were actually built, but for the 
objectives cf this analysis they are left as continuous real 
variables. for example, the number of shells may have an 
integer value in the actual system but this requirement can 
be accomodated by solving the system for a real value then 
integerizing the number of shells and solve the problem with 
the next higher and the next lower integer number of shells. 
The lower cf the two solutions would be chosen. The final 
rescluticn of the correct number of shells would 
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Table 4 

DEilKIIICN OF THE ELEMENTS OF THE VECTOR X 



X (1) 

X (2) 

1 (3) 
x(4) 
x(5) 

X (6) 

x(7) 

x(8) 

x(9) 

x(1C) 

1 ( 11 ) 

1 ( 12 ) 

1(13) 



velocity of seawater in the boiler 
” ” '• '• “ condenser 

cutlet pressure in the boiler 
inlet '• '• ” condenser 

tube length in the boiler 

" ” ” ” condenser 

aass flow rate of the working fluid 
diaaeter of the boiler shell 
“ '• ” condenser shell 

outside dia. of a tube in the boiler 

II II II II If II II condenser 

number of boiler shells 
'• ” condenser shells 
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probably depeod on other than thenaodynaaic and cost 
considerations, such as a reguirement that cne extra shell 
be included to allow for one shell being down for 
maintenance. Obviously, the number of tubes in a shell must 
be an integer value but one extra tube cut of several 
thousand is a minor matter. In addition, there may be a 
requirement for extra tubes to allow for a percentage to be 
plugged when they develop leaks so that the entire shell is 
net put cut cf commission by a minor leak in one tube. 

1 • Objective fu nct ion and the C onstra ints 

Ihe objective function is the sum of the costs of 
the boiler and the condenser, and the pumps. The model 
assumes that there is one circulating pump and one feed pump 
for each shell in the boiler and that there is one 
circulating pump for each shell in the condenser. Since 
there is no cost analysis for the turbine, no assumption is 
needed about the number of turbines in the system. 

There are nineteen constraints in the present 
pregram, cne cf which is optional. The first three ccncern 
the saturation pressures in the boiler and the condenser. 
The next twelve place restrictions on various dimensions of 
the boiler and the condenser. The last three require that 
seme intermediate calculations stay positive. 

Using the objective function as a constraint has 
proved advantageous in cases where the problem was 
ncn-ccnvex or where the program had difficulty in finding a 
solution. It is used by setting a boundary on the objective 
function so that the constraint is either infeasible or 
slightly feasible at the starting point. The proper choice 
must be round by trial and error. One caution to note is 
when equality constraints are involved. If the starting 
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pcint is neat ths optimum, the program tries to move away 
from the boundary and the eguality constraints may not be 
satisfied wten the program finishes. 



The need for the ether eighteen inequality 
constraints is not obvious when the problem is first 
formulated. The thermodynamics of the problem require, for 
example, that the propane saturation temperature be less 
than seawater temperature in order for the bciler to absorb 
energy, bet the program has no way of knowing that, unless 
it it specified. If the program can improve the cost and 
satisfy the ether constraints by reversing the relaticnship 
it will dc so. This inversion may take place only for a 
short time during the solution of the problem, but, if the 
equations fail under this condition, the censtraint cn the 
temperatures in the bciler is necessary. Most of the other 
constraints are needed for the same reason. 

Seme constraints are not added to the problem since 
they are net binding in the solution. There are twe such 
constraints. The first is the segmental area required to 
allow separation of the liquid droplets from the vaper. The 
segmental area set aside by the 40X rule-of-thumb is more 
thar that required by the consideration of the physical 
properties cf prepane. The ether censtraint is the 
restriction cn the maximum heat flux by Mostinski's 
correlation, equation (37) . The constraints dc not become 
binding in the solution. 

The three eguality constraints contain the 
calculations cf the properties of the various components 
modeled in this analysis. The first cne concerns the 
bciler. It states that the actual energy absorption rate of 
the boiler shells must be equal to the energy absorption 
rate required by the thermodynamic cycle. Similarly, the 
second eguality censtraint states that the heat rejected by 
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th€ condenser must egual that required by the thermodynamic 
cycle for the condensation of the vapor from state point 5 
back to state pcint 1 as shown in Fig, 3, The third 
constraint calculates the power required by the various 
pumps and the pcwer extracted from the turbine to run the 
system's pumps and to produce electrical power for 
consumption. This constraint requires that the power cut of 
the turtine equal the power used by the pumps and the 
generator. The net useful power cut of the system is a 
design objective and is specified as 25MW herein, 

2, The Cp tim ization Method 

Since the arrival of the electronic computer, a 
great many nonlinear programming techniques have been 
developed frcm early theories and more current research. 
The computer program, called the Sequential Dnconstrained 
Minimization Technique (SOMT) , developed by Hylander, 
Helens, and McCormick in Bef. 22 is used to optimize the 
model. Their program implemented much of the theory 
contained in Ref. 23, by Fiacco and McCormick, about 
nonlinear programming using unconstrained minimization 
techniques, SDMT-Version 4 is chosen because of its ability 
to solve a wide variety of problems. The thecretical 
requirements for SOMT to be guaranteed to find a local 
minimum is that there exists some point which satisfies the 
inequality constraints and that there not be a local minimum 
at points where X goes to infinity. If certain convexity 
conditions are satisfied, then the local minimum found is a 
global minimum. 

£CMT uses a penalty function tc transform the 
constraired cptimization problem 
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oicimize 



f (X) 

subject tc: g (X) > 0, i=1,2, 

i 

h (X) = 0, i=ai+1 , . . . . ,m+n 
i 



into a sequence of unconstrained problems of the form 



m m+n 

P(X,r) = f<X) - r Z In g (X) ♦ 1/r I (h (X) ) 2 (61) 

x=1 i i=m+1 i 

where E(X,r) becomes the function to be minimized. The 
inequality constraints, g (X) , cause a large penalty to cccur 
as a boundary is approached. The equality constraints add a 
large penalty whenever h (X) departs from zero in either 
direction. The size of the penalties due tc g (X) and h (X) 
is contrclled by the arbitrary pa-rameter "r". The algorithm 
solves the problem by solving a series of subproblems. For 
r^ >r^>rg . . .>0, SDMT minimizes P(X(r„),r,^) in each subproblem 
in order to find the minimum of f (X) as "r” tends to zero, 
ahen a minimum of a subproblem is found, "r" is reduced by a 
factor assigned by the user and the new subprcblem is 
solved. The process continues until the convergence 

criteria chosen by the user is satisfied. Figure 7 
illustrates schematically the results of the sequence of 
suhprcbleits and the qualitative shape of the . penalty 
function. Three methods of minimization are included in the 
program. The user may choose either the generalized 

Newton-Eaphscn, steepest descent, or McCormick modification 
of the Fletcher-Powell method. 



There are several other options that control the 
operation of SOMT. The user may cheese the initial 
penalty, r, cr, if the problem saxisfies certain conditions, 
the user may elect to have SOMT compute the initial "r'*. 
There are three choices for the convergence criterion cf the 
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FIGURE 7. PROBLEM SOLUTION SEQUENCE 
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pECblem. Ihe criterion used depends on the inequality 
penalty ten being close enough to zero. Another option is 
the choice of the calculation of the subproblem stepping 
criterion . Two of the choices compute the inverse cf the 
secend partial matrix and the third choice computes only the 
gradient of the penalty function (Vj^P(X,r)). One cf the 
Bcst useful options is extrapolation. After two or three 
suhproblens have been solved SDMT is able tc make a first or 
secend creer extrapolation to the starting point for the 
next sukprcblem. This can speed up the solution cf the 
problem measureably. 

The characteristics of the problem limit the choices 
for each cpticn. Ihe penalty factor couldn't be set by SOHT 
due to the presence of equality constraints and because the 
initial 1 is usually not close to any boundary. The stepping 
criterion fer the problem is that the inequality penalty 
term become less than some value, usually 10“'^. 
Extrapolation is used, but it makes little difference in the 
answer cr tc the solution time whether a twe point or three 
point estimate is made. The subproblem convergence 
criterion is the magnitude of the gradient of E(X,r) 
becoming less than some small number. 

Ihe choice cf starting penalty, r, affects 
significantly the solution time and even whether a solution 
is found it seme problems according to Bef. 7. If the 
initial "r'* is toe small, the subproblem optimum may be the 
problem's optimum which may be hard to find cr require more 
calculations than if a much larger '*r“ is chosen. If "r" is 
toe large, the first few optima lie near the center of the 
feasible region and are unrelated to the optimum cf the 
objective function. 

Scaling of the variables, constraints, and the 
objective function is necessary. The variables must be 
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witfcin tKC cr three orders of magnitude of each other. For 
example, the working fluid flow rate is of the order of 10^ 
whereas the tube diameters are of the order lOo. If the 
finite difference interval is 10-^, the partial of the 
eguations with respect to the mass flow rate may be lost due 
tc the roundoff error providing the only significant 
figures. Rith the constraints and the objective function, 
unsealing may result in one term in the penalty function 
being dcainant, and the algorithm finds that most cf the 
decrease in the penalty function takes place in the dominate 
term. Ihe result is that SOMT works hard to decrease that 
tern and does not see the other terms resulting in a 
meaningless solution. 

In making first and second order estimates of the 
solution, SCMT may generate negative numbers in the X 
vector. Hhen it evaluates the constraints and objective 
function for -this point, the equations fail and may cause 
the program to fail depending on the severity of the failure 
and how many times the error has cccured. This problem is 
circumvented by testing the vector upon entry into the 
subroutine where the objective function and the constraints 
are evaluated and setting the constraint being evaluated 
equal to seme negative constant. 

Ihe main program of SUMT is altered so that the 
problem can be restarted at the current solution without 
reading in a new problem. This is necessary because the 
pregram would give all the indications of having feund a 
solution, hut if the problem is restared it could often find 
a lew value for the objective function. The new solution is 
compared to the previous solution and if the two solutions 
are not close epough the problem is restarted after making 
seme changes to particular parameters, options, or 
constraints depending on the objective the user has in mind. 
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III. RESDLIS AND CONCLO SIGNS 



One cf tbe majcr tasks in the development of a nonlinear 
model is the linking of the model with the optimization 
program in order to analyze the model. The selection cf the 
proper nonlinear programming technique is essential tc the 
solution. Several techniques require that the objective 
function and the constraints be in a particular form. For 
example. Geometric Programming developed by Cuffin, Peterson 
and Zener, in Ref. 24, requires that the equations be 
posyncmials which is not true of the model developed in this 
paper. SOMI was chosen because cf its wide applications to 
problems cf anj type. It only requires that there exist 
some feasible point and that local minima not exist at 
infinity. It has successfully solved many ncnconvex 
problems. 

Twc versions of the analytical model are presented. The 
first is called the iterative version because the 
calculation cf the heat absorption rate and heat rejection 
rate require iterating on some quantity to find them. 
Hostinski's correlation, equation (36), contains the heat 
flux which is an unknown. It is not possible to solve the 
system cf equations explicitly for the heat flux; therefore, 
a successive approximation technique is used. The 
calculation cf the heat rejected in the condenser involves 
the same process with the temperature of the outside tube 
wall as the iterating variable. SDMT is urable to arrive at 
a consistent solution in the time allowed, up to 90 minutes 
of computer time. The solutions are quite different for 
each try with nc pattern • evident about tbe direction the 
scluticn might lie. 

✓ 
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There are several possible explanations for this 
behavior. Perhaps the behavior of this version of the 
problem is such that SOMT is not the best choice fcr the 
optimizinc algorithm or possibly the problem cannot be 
solved at all in this form. Another source of problems 
could be the numerical differentiation procedure. Numerical 
differentiation is subject to severe truncation and 
round-off error. The subproblem always terminates when SUMT 
is unable to reduce the penalty function during the next 
iteration, from one estimate of X to the next estimate. 
Remember the convergence criterion, never satisfied for this 
version of the problem, is the magnitude of Vj^E(X,r) 
becoming less than some small number. No conclusive results 
could be obtained with this version so the second version is 
used fcr the present analysis. 

In this second version the heat transfer coefficients on 
the Propane side of the tube are fixed. This enables the 
eguality constraints for the boiler and the condenser tc be 
calculated to the full accuracy of the computer without 
increasing the time reguired. The correlations, eguation 
(36) and (48) , are used for comparison with the fixed value. 
If the calculated value and the fixed value are within 20 to 
2531 of one another, the solution would be considered to be 
good since this is the tolerance claimed for the 
ccrrelaticns . 

SDMT seemed to be able to solve the new prcblera because 
a series of seven solutions have been produced where the 
values of the elements of X and of the objective function 
vary nc more than 3.735 and ^. 2 %, respectively, about the 
mean when different starting points and different SUMT 
control parameters are selected. See Table 6 for a summary 
of the scluticns obtained from the assumed data, as shewn in 
Table 5. The design implications of these results are quite 
reasonable. The range shows the spread in the solutions. 
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Table 5 

PROBLEM INPOT DATA 



Eouliag Thermal Resistance 
Eciler and Condenser 
Elevation cf Bciler 
Atcve Condenser 
Channel Allowance 

Efficiencies for the Turbine 
and the Pumps 

Seawater Temperatures - Hot 

- Cold 



Net Output Power 

Eundle factor (F) 

Cost Icdexes (I) 

Tube Hall Thickness 

Min. Cistance Between Tubes 

Heat Transfer Coefficients 
fcr Eciler and Condenser 



0.0 hr-op/BIO 
(0.0 °C/kH) 

15 ft 
(4.57 m) 

0.83 ft 
(0.25 m) 

8551 

75op 

(23.9°C) 

40°F 

(4.4<>C) 

8.532x107 BTO/hr 
(25 MH) 

3.0 

1 .00 fcr all ccmp . 
0.028 in 

(0.071 cm) 

0.5 in 
(1.27 cm) 

145 BID/ft2-hr-0F 
(0.822 kH/m2-oc) 



Materials 



Boiler and Condenser 

Feed Pumps 
Circulation Pumps 



Carbon Steel Shells 
Brass Tubes 
Bronze 

Stainless Steel 
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Table 6 



Variable 


EXAMPLE PSOBLSa 
Banqe 


Seller : 


Velocity 


3.323-3.347 


Eressure 


(1.013-1.020) 
aiprox 0 


lute length 


39.16-40 .00 


Shell dia. 


(11.91-12.19) 

24.12-25.00 


lube dia. 


(7. 352-7.620) 
0. 756-0.761 


# of shells 


(1. 919-1.934) 
7. 141-7.409 


Cerdenser; 


Velocity 


2.793-2.823 


Eressure 


(0.851-0. 860) 
91.53-91.60 


lube length 


(63 1. 1-631.6) 
39.35-40 .00 


Shell dia. 


(11.99-12.19) 
44 .69-45 ,00 


lube dia. 


(13.62-13.72) 

0.7632-0.7686 


1 of shells 


(1.938-1.952) 

3.297-3.401 


Mass flow rate 


22. 12-22.21 


System cest 


(10.03-10.07) 

7.578-7.673 



3OL0TIOH SUHflABY 



Mean 


Var 


Units 




-5- 




3.33 5 


0.7 


ft/sec 


(1.016) 




(a/sec) 


118.2 


0.0 


lb/in2 


(815.0) 




(’<P) 


39.71 


2.1 


ft 


(12.10) 




(m) 


24.91 


3.5 


ft 


(7.593) 




(m) 


0.758 


0.7 


in 


(1.926) 




(cm) 


7.228 


3.7 


— 


2.808 


1 . 1 


ft/sec 


(0.856) 




(ra/sec) 


91.56 


0.1 


lb/in2 


(631.3) 




(JtP) 


39.78 


1.6 


ft 


(12.12) 




(m) 


44.89 


0.7 


ft 


(13.68) 




(m) 


0.7660 


0.7 


in 


(1.946) 




(cm) 


3.331 


3. 1 


- 


22.17 


0.4 


10 ®lbra/hr 


(10.06) 




(10*kg/hr) 


7.607 


1.2 


106 $ 
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The arithaetic meaQ of the seven solutions is calculated. 
The difference betueen the naximum and the ainimum values is 
shewn as a percent of the aean. The percent variation gives 
scae idea of how well the solution was determined. Note 
that the tube length and shell diameter constraints are 
binding for both the boiler and the condenser. The 
subproblem still terminates only when the penalty function 
does not decrease from one iteration to the next. The 
magritude of ^P(X,r) is not becoming small enough to 
satisfy the convergence criterion. The value never becomes 
less than C.02 and is usually greater than 0.5 when a 
subproblem is stopped by the program. 

Useful results are obtainable, however. The engineer 
can learn much about his problem even from a single point. 
For example, from the cost breakdown in Figure €, the 
portion of the costs attributable to each component can be 
seen. The beat exchangers, of course, dominate with 783? of 
the cost. The circulation pumps are about 21 % of the cost 
with the feed pumps being T.235 of the total cost. Since 
the efficiency assumed for the circulation pumps is prebahly 
toe high, these pumps would be even mere costly in terms of 
energy and mcney. This could indicate the priorities that 
should be set on the design efforts for improvements cf the 
various cempenents with the highest priority being the 
highest percentage cost component. 

The system’s power distribution is shown in Table 7. 
The sigrificant features are the large energy absorption 
rate for the small useful output power and the parasitic 
pumping power. The lower the efficiency the greater must be 
the power absorbed for the same output. This is to be 
expected for this kind of system. The parasitic pumping 
power (1231 of the turbine output) is about twice what it is 
for a fessile fuel plant. 
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CONDENSER 

$3,333,000 

43.91 



BOILER 

$2,528,000 

33.3% 




FEED / BOILER $978,000 

PUMPS L CIRCULATION 12.9% 

$94,100 PUMPS 

1.2% $654,000 



8 . 6 % 



FIGURE 8. COST BREAKDOWN 
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Table 7 

SOIDTICN POHEE BOCGET 



Beat absorption 
Turbine Shaft Power 
Eeed Pumps 

Elr Circulation Pumps 
Ccnd Circulation Pumps 
Eeat Bejected 
Ket Cseful Power 



BTO/hr 


Matts 


3.56x109 


1042.8 


9.69x107’ 


28.4 


4.49x106 


1.3 


2.58x106 


00 

t 

o 


4.51x106 


1 .3 


3.47x109 


1015.7 


8.53x107 


25.0 
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Cdc€ the model has been optimized with the number of 
shells as real variables, the problem oan be re-solved twice 
by fixing the shells at the next higher and the next lower 
integer numfcer. For the higher number of shells, the cost 
increased slightly as expected, from 7.607 million dollars 
average to 7.983 million dollars. The problem with the 
lower integer number of shells could not be solved due to 
the constraints cn the tube length and the diameter cf the 
shells being binding before the shell numbers were 
integerized. The program is unable to find a new solution 
that satisfies the equality constraints. 

Table 8 shows some miscellaneous data from the model 
that deserves comment. The overall heat transfer 
coefficient cf the boiler is approximately three times that 
for the condenser which is the major reason for the lower 
cost of the boiler relative to the condenser. However, 
remember that the bundle factor, equal to 3, is a part of 
the heat transfer coefficient. More data must be gotten 
abcut this factor in order to improve the reliance on any 
solution. There is a significant difference between the 
Prepane heat transfer coefficients that were assumed ard the 
theoretical values calculated from the correlations. Other 
values should be tried until the fixed value is withir the 
2051 tolerance on the theoretical value. The flew of 
seawater through the tubes in both the boiler and the 
cerdenser is turbulent. Apparently, the resulting cost 
savings due to the improved heat transfer rate through a 
smaller area is more than rhe increased costs of the higher 
pumping power needed. It seems likely that the solution 
obtained is globally optimal; several different starting 
points were tried. The problem has been started in both the 
laminar and the turbulent flow regions with no change in the 
results. 
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Table 8 

ADDITIONAL DATA PROM SAMEIE 
PROBLEM SOLUTION 

Condens er Uni ts 

Scrking fluid Heat Transfer Coefficients 



Assumed 


145. 


145. 


BTD/hr-f t2-op 




(0.822) 


(0.822) 


(kW/m2-cC) 


Theoretical 


227. 


241. 






(1.29) 


(1.37) 




Surface Area per Shell 


2.39x105 


1.27x105 


ft2 




(2.22x10+) 


(1.18x105) 


(m2) 


Seawater Cutlet Temp 


72.3 


43.8 


op 




(22.4) 


(6.55) 


(°C) 


Shell Saturation Temp 


66.7 


49.9 


op 




(19.3) 


(9.94) 


(°C) 


Seawater Mass Plow Rate 


4.39x10« 


9. 11x10« 


Ifam/hr 




(1 .99x10«) 


(4. 13x109) 


(kg/hr) 


Mean Temp Difference 


6.67 


7.85 


Op 




(3,71) 


(4.36) 


(°C) 


Overall Beat Transfer 


310.6 


105.0 


ETU/hr-ft2-op 


Coefficient 


(1.76) 


(0.595) 


(kH/m2-oc) 


Heat Plui 


2070.0 


824.9 


ETU/ft2-hr 




(6.528) 


(2.602) 


(kW/m2) 


Seawater Bejrclds to. 


19,800. 


9,960. 


- 
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IV. R EC OH HEMATIC NS 



The Ei^€£€nt model has a great deal mere infomaticn to 
offer, but seme improvements in the model shculd be made to 
improve the ability of a nonlinear program to solve the 
problem. The first recommendation would be to study the 
equality constraints and place inequality constraints on any 
variable or grouping of variables if they both can become 
negative and also must have their logarithms computed. 
Examples, already included in the program, are the number of 
tubes in the shells. The heat transfer coefficients for the 
working fluid side in the boiler and the condenser could be 
included as equality constraints. This can be dene by 
adding the heat transfer coefficients to the vector X and 
requiring that the variables be equal to the coefficients as 
computed by the appropriate correlation. The same could be 
done with Ealen and Small's property-dependent calculation 
of the separation area by equation (12) and (13). From an 
overall standpoint, the results seem to indicate either the 
problem is too difficult to be solved efficiently it its 
present forv or that another nonlinear program, other than 
SDHI, might he better suited to the problem. Numerical 
differentiation of the constraints in this problem do not 
seem to help the program to find the solution since the 
magnitude of ^P(X,r) is not near zero when a subproblem 
steps. 

Once the optimization scheme is functioning properly, 
the model can be used to study the effects of fouling and 
corrosion and the costs of their prevention. Many proposals 
have been made that must be analyzed in the context of the 
total sjsten. Seme such as Amertap, where plastic foam 
balls are pumped through the tubes to clean the heat 
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transfer surfaces, add a cost penalty frcn the capital and 
cperating ccst of the ftmertap system and from the increased 
pumping pcuer required. Coatings add the capital costs of 
the coating and the added capital cost of increased heat 
transfer area required because of the increased thermal 
resistance cf the coating. 

A sensitivity analysis can be performed on the model to 
test the assumptions and the various ccnstants. The 
sensitivity to changes in these factors can indicate how 
much one should be willing to pay for the improvement 
sought. High sensitivity indicates areas where the greatest 
effort for improvement should be desired or the priorities 
should he set cn research into the unknown areas cf the 
problem. The trade-off between cost and performance can be 
studied by changing the output power and looking fcr the 
lowest ccst per kilowatt-hour, which is the measure by which 
ccnventicnal power plants are compared. One could look into 
the cost and the benefits of increasing the depth from which 
the cold water is pumped versus the improved temperature 
differential possible from using the colder water. 

For the capital cost of the system to be complete, the 
turbine design must be included and ccst estimates made. 
The turbine will be unique due to the very low pressure 
potential drop it must use. Some work has been done by 
other groups working on OTSC, Hef. 25 and 26. fieference 27, 
by ficbert I. Bartlett, can provide costing and design 
infcrmaticn although he deals only with steam turbines. The 
remaining part of the system that is desireable for 
inclusion it the analysis presented herein is the cold 
seawater intake pipe. It will be on the order of 1,0C0 ft 
(3C4.6m) Iccg and 100 ft (30.4 m) in diameter. 
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